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I. Introduction 



In recent years cosmological models with a scalar field acquired the great popu- 
larity because they can serve as a most natural basis for the inflationary cosmology 
1 . Indeed, the presence of scalar field in the model under consideration provides 
us with the existence of effective cosmological constant at the early stage of the 
cosmological evolution and with an opportunity to describe the decay of this con- 
stant and the "graceful exit" from inflation and transition to the Friedmann stage 
of evolution at a proper moment. 

On the other hand scalar field is an integral part of modern models in parti- 
cle physics. Moreover, the main part of papers devoted to quantum-cosmological 
description of the quantum origin of the Universe and to the definition and con- 
struction of the wave function of the Universe consider the models including scalar 
field which after the "birth" of the Universe is driving inflation 2 ~ 4 . 

Side by side with the comparatively simple models based on the simple La- 
grangians of the scalar field were developed rather complicated schemes considering 
the non- minimally coupled scalar field 5 , complex scalar field 6,7 or the scalar field 
combining complexity and non-minimal coupling 8 ' 9 . 

However, even the dynamics of the simple cosmological model, including gravity 
and minimally coupled scalar field with simple potential including only massive 
term is rather rich and deserves studying. The dynamics of the minisuperspace 
cosmological models with the massive real scalar field for the flat, open and closed 
Friedmann universes was studied in papers 10 ' n in terms of phase space and theory 
of dynamical systems. It was noticed that the dynamics of closed model (which is the 
most interesting from the quantum-cosmological point of view) is more complicated 
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than that dynamics of open and flat models. This dynamics allows the transitions 
from expansion to contraction and the existence of points of maximal expansion and 
minimal contraction in contrast with cases of open and flat cosmologies. Moreover, 
closed spherically symmetric models cannot expand infinitely and should have the 
points of maximal expansion provided the matter in the model under consideration 
satisfies the condition of oncrgodominance 12 . The presence of points of maximal 
expansion and minimal contraction open the possibility for the existence of the 
trajectories of evolution of the Universe escaping singularity and oscillating in a 
periodical 13 or in an aperiodical 14 way between turning points. The possibility of 
existence of such trajectories or, in other words, non-singular universes filled with 
scalar field was discussed also earlier in Rcf. 15. 

Here we would like to consider the simplest cosmological model with massive 
real minimally coupled scalar field without self-interaction. We shall consider closed 
Friedmann model in minisuperspace including only two variables - cosmological ra- 
dius a and homogeneous mode of inflaton scalar field <p. Because all the trajectories 
in such a model have the point of maximal expansion one can give the classifica- 
tion of the trajectories starting their evolution since these points. Studying such 
a classification we ignore the "prehistories" of the trajectories under consideration, 
i.e. we do not study their evolution before acchieving points of maximal expansion. 
Such an approach simplifies the classification, because it gives us an opportunity 
to treat in the same way all the different kinds of trajectories - trajectories which 
exist during some finite intervals of time, i.e the trajectories which were born in 
singularity and disappear in it, the infinitely oscillating trajectories escaping falling 
into singularity and two kinds of "semi-infinite" trajectories which were either born 
in the singularity or disappear in it. 
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Suggested classification according to points of maximal expansion is convenient 
because localization of these points (i.e. a set of possible points of maximal expan- 
sion as well as those of minimal contraction) is well known and was described in Ref. 
10 in terms of phase space and in Ref. 14 in terms of configuration space. General 
formulas describing such points for more wide class of models were presented in Ref. 
9. 

We shall distinguish trajectories which are monotonical in a and ip i.e. tra- 
jectories which fall to singularity with monotonically changing value of scalar field 
and the trajectories possessing bounces i.e points of minimal contraction when a 
is equal to zero and trajectories having extrema in the value of scalar field ip - 
we shall call them y>-turns. It is remarkable that there is some regularity in the 
localization of the points of maximal expansion corresponding to different kinds of 
trajectories. The points is that the regions corresponding to trajectories falling to 
singularity after some definite number of </?-turns are separated by regions corre- 
sponding to bouncing trajectories. This phenomenon was detected by us due to 
numerical simulations and will be presented in detail in the second section. 

Third section will be devoted to comparison of our model with even more simple 
model with massless scalar field 16 . This model is of interest because it is exactly 
solvable and have only trajectories without bounces and tp-turns. So it is interest- 
ing to try to design some kind of perturbation theory for our model with massive 
scalar field where exactly solvable equations of motion for massless model play role 
of zero-order approximation. Such scheme seems promising because of interest to 
exactly solvable cosmological models rising in recent years 17 . We shall see that for 
our model one can design such perturbative scheme and even to get the first-order 
corrections in an explicit way. However, this perturbation theory does not work in 
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the vicinity of singularity. Nevertheless, it will be shown that combining the solu- 
tions of equations of motion obtained in the framework of perturbation theory with 
some general properties of solutions of full (massive) theory one can obtain some 
useful information concerning different regimes of cosmological evolution. Perhaps, 
the approach developed in the present paper can be useful for analysis of the dy- 
namical behavior of cosmological models which are close to some exactly solvable 
cases. 

II. The model, equations of motion and their numerical investigation 

We shall consider the cosmological model with an action 

■,2 



S = J <£xj—g | ^-R + \<rdp<pd v <p - imV 



(2.1) 



For the closed Fricdmann model with the metric 

da 2 = N 2 (t)dt 2 - a 2 (t)d 2 n (3 \ (2.2) 

where a(t) is a cosmological radius, N - a lapse function and d 2 fl^ is the metric 
of a unit 3-sphere and with homogeneous scalar field ip the action (2.1) takes the 
form 

_ I" „ I Jim 2 , / n \ 2 1 tn^ m^in^ 1 

(2.3) 



S = 2tt 2 J dtNa 3 I ^ 



8tt 



\Na) + a 2 



ip 2 m 2 (p 2 
2N 2 2~ 



Now choosing the gauge N — 1 we can get the following equations of motion 



• 2 i \ -2 2 2 

l p ■■ a 1 \ ap> m ip a 



'd-\ 1 )+— 1 — =0, (2.4) 

16tt V 2a 2a J 8 8 ' V ; 

i/5+ — +m 2 (p = 0. (2.5) 

Besides, we can write down the first integral of motion of our system 

3 a 2 
~—m 2 P (d 2 + l) + — ((p 2 +m 2 p 2 )=0. (2.6) 

S7T Z 



It is easy to see from Eq. (2.6) that the points of maximal expansion and those 
of minimal contraction, i.e. the points, where d = can exist only in the region 
where 

2 3 m 2 P 

which represents the field in the half-plane < a < +oo, — oo < ip < +oo restricted 
by hyperbolic curves ip < \J^ 1 ^ an d <p > ( see Fig. 1). Sometimes, 

region defined by nonequalitics (2.7) is called Euclidean or "classically forbidden". 
One can argue about validity of such definition (see, for detail Ref. 9), but we 
shall use it for brevity. Now we would like to distinguish between points of minimal 
contraction where a — 0, a > and those of maximal expansion where a = 0, a > 0. 
Let us put a = 0, in this case one can express <p 2 from Eq. (2.6) as 

^ = £^-»V. (2.8) 

Substituting (2.8) and a = into Eq. (2.4) we have 

4irm 2 p 2 a 2 



m. 



(2.9) 



^From Eq. (2.9) one can easily see that the possible points of maximal expansion 
are localized inside the region 

^<^^V (2-10) 

while the possible points of minimal contraction lie outside this region (2.10) being 
at the same time inside the Euclidean region (2.7) (see Fig. 1). 

We shall consider the trajectories beginning at some points of maximal ex- 
pansion in the region defined by (2.10). As it has been already mentioned in the 
Introduction we shall trace out only the parts of the trajctories beginning since the 
points of maximal expansion ignoring their "prchistories" . 
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For the velocities ip defined by Eq. (2.8) we choose direction "up", i.e. positive 
sign (the picture with velocities ip pointed out "down" can be obtained by mirror 
reflection in respect to axis <p = 0). Numerical investigation shows that for the 
points of maximal expansion placed in the region denoted in Fig. 2 by number 
trajectories monotonically approach to singularity a = 0, ip = +oo (some typical 
trajectories of this kind are shown in Fig. 3a). It looks quite natural because in the 
region close to the ordinate axis a = behavior of our system is close to that for 
the system with massless scalar field where bounces and ip - turns are absent. Then 
trajectories beginning in the region denoted in Fig. 2 by number 1 have a bounce. 
The typical trajectories of this kind are shown in Fig. 3b. Then we have region 
which we shall denote by number 1. In this region we also have a bounce, but before 
this bounce we have ip - turn. The typical trajectories beginning in the region 1 are 
shown on Fig. 3d. The boundary between the regions 1 and 1 is the place where 
some periodical trajectories escaping singularity can have their points of maximal 
expansion (such trajectories were first mentioned in Ref. 13 and analyzed in detail 
in Ref. 14). In particular, through the point ip = lying on the boundary between 
the regions 1 and 1 goes the symmetric in <p periodic trajectory depicted in Fig. 3c. 

In the region 1' are localized the trajectories which have not a bounce, but ip - 
turn after that they tend to singularity a = 0, ip = — oo (see Fig. 3e). The boundary 
between the regions 1 and 1' corresponds to trajectories where the points of bounce 
(a = 0, a > 0) degenerates into inflection points (d = 0, a = 0). 

Then in the region 2 (see Fig. 2) we have the points of the maximal expansion 
for trajectories which after going through t^-turn have a bounce at ip < escaping 
in such a way the fall into singularity (see Fig. 3f). In the region 2 we have the 
trajectories which have a bounce in the lower half-plane, just after two t^-turns. The 
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boundary between the regions 2 and 2 again contain trajectories including turning 
points where we have simultaneously d = 0, (p = 0. Some of these trajectories are 
periodical. One of such trajectories is shown in Fig. 3g. 

Then starting since the point of maximal expansion in the region 2' (see Fig. 
2) we have the trajectories without bounces, but with two <p-turns (one at ip > 
and one at tp < ) after which they go to singularity a — 0, tp — +oo (see Fig. 3h). 

Now one can easily understand what kind of trajectories corresponds to the 
points of maximal expansion localized in regions 3, 3, 3', 4, 4, 4' and so on. 

It is worth noticing that the most peculiar form has the boundary between 
the regions 1 and other regions in the upper half-plane (ip > 0). This boundary 
consists from two curves (see Fig. 2). The left one is finite and divide the region 1 
from the region 1 while the right curve is an infinite one and goes almost parallclly 
to the hyperbolic curve separating the possible points of maximal expansion from 
the possible points of minimal contraction (see Fig. 2). It is interesting that this 
upper curve touches all the regions beginning since T. 

In Fig. 3i is shown the trajectory whose point of maximal expansion is placed 
in the upper branch of the region 1. This trajectory immediately after the bounce 
has ip - turn, then goes through the second point of maximal expansion falls into 
singularity in the lower half-plane. 

All the regions beginning since T are compact in contrast with regions and 1. 

Now one can pay attention to a special class of trajectories having the points of 
"full stop" , i.e. the points where a — 0, tp — simulataneously. These points of full 
stop can lie only on the boundary between Euclidean and Lorentzian region. Nat- 
urally, the trajectories beginning with zero velocities at this boundary constitute 
the tiny minority in the set of all possible trajectories. Nevertheless it is interesting 
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to trace out the correspondence between the points of full stops of these trajec- 
tories and their points of maximal expansion. Here it is necessary to remember 
that in every point of maximal expansion we have two possible directions for <p. 
Our classification was constructed for the direction "up" . The structure of regions 
corresponding to the direction "down" can be obtained by reflection of picture in 
respect to axis ip = 0. We shall distinguish between regions "up" and "down" by 
symbols f and J. correspondingly. 

Thus, let us consider the trajectories which begins with zero velocities a = 
0, p = in the upper half-plane ip > 0. As was explained earlier the corresponding 
trajectories can have their points of maximal expansion on the curves separating 
the regions k from the regions k. Besides these boundaries the points of maximal 
expansion can be disposed on the infinite curve separating the region 1 from all other 
regions. Indeed, numerical investigation shows that the trajectories beginning at 
< ip < ipo , go into Euclidean region and have the first point of maximal expansion 
just on this infinite curve separating the region If from other regions Here (po is the 
point where the direction of the motion at the initial moment coincides with the 
direction of the tangent to the curve given by an equality in (2.7) , i.e. the point 
where 

l = d f. (2.11) 

a da 

Using Eqs. (2.4), (2.5),and (2.7) one can find 



Then we consider trajectories beginning with zero velocities at <p > ipo. Moving 
up the beginning of trajectory along the hyperbolic curve separating Euclidean and 
Lorentzian regions beginning since po we shall have the corresponding point of 




(2.12) 
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maximal expansion moving down along the boundary between the regions If and 

1 j. Then, after approaching the final point of this curve (i.e. the crossing between 
this boundary curve and the curve separating the points of maximal expansion 
from points of minimal contraction) the closest to the beginning of motion point of 
maximal expansion jumps to some point of the boundary between regions 2J. and 

2 J.. Then while the beginning of the trajectory is moving up along the hyperbolic 
curve separating Euclidean and Lorentzian regions the corresponding point of max- 
imal expansion is moving up to until the moment when it approaches to the curve 
separating the points of minimal contraction and maximal expansion in an upper 
half-plane. Then it jumps to some point on the curve separating the regions 3| 
and 3 f and so on and so forth. Thus we have seen that every point on the upper 
branch of the hyperbolic curve separating Euclidean and Lorentzian regions has a 
counterpart on one of curves separating regions k and k and the intervals denoted 
with sign "f" is interchanged with intervals denoted by that of "J." . The picture for 
the down part of hyperbolic curve disposed at ip < will be quite the same with an 
only difference that signs f should be substituted by ones J. and vice versa. 

Now let us turn back to the analysis of the general set of trajectories in our 
problem. The most interesting for us will be trajectories which have bounces, and 
especially the trajectories which have a lot of bounces and as an ideal situation an 
infinite number of bounces (i.e. the trajectories escaping singularity). 

Let us look at the trajectories whose points of maximal expansion are localized 
in the regions 1, 2, 3 . . . , or in the regions 1,2,3,... i.e. the trajectories having at 
least one bounce. We shall study the structure of these regions from the point of view 
of localization of their second point of maximal expansion. The substructure of the 
region 1| is represented on Fig. 4, while the forms of the trajectories corresponding 
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to various subregion of regions 1 and 1 are presented in Fig. 5. 

On the right near the boundary with the region 1 we have subregion If 1 J. 
which corresponds to the trajectories which have after a bounce the second point 
of maximal expansion in the region 1 j after which they have the second bounce. 
The next subregion to the left from the subregion If 1J. is the subregion If l'j 
which corresponds to the trajectories which after bounce have the point of maximal 
expansion in the region l'j. After this point of maximal expansion they have a ip - 
turn and then fall to singularity a = 0, ip = — oo . 

Then we have the subregion If 1 j which corresponds to trajectories with two 
bounces. Then we have subregion If l'f . Then follows the subregions If 2f , If 
2f, If 2'f , If If, If 2'f , If 3f and so on and so forth. 

The structure of region If is much more simple: we have here only two subre- 
gions If If and T Of. 

Thus inside the regions If, If we have an infinite set of subregions whose 
arrangement approximately repeat the structure of the whole field of the possible 
points of maximal expansion presented in Fig. 2. 

One can see that the structure of couples of regions 22; 3,3; 4,4 ... is quite 
similar to that of couple of regions 1,1. 

Studying the structure of subregions corresponding to two bounces one can 
find inside each of them the same "subsubstructure" of infinite set "subsubregions" . 
Continuing this procedure ad infinitum one can observe that the field of the local- 
ization of the points of maximal expansion corresponding to trajectories escaping 
singularity can be find as a result of an infinite procedure at each stage of each we 
encounter self-similar structures. Such a self-similarity of structures appearing at 
different scales points out on the fractal nature of the set obtained as the result of 
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infinite procedure 18 . Thus, while the set of trajectories escaping singularity and 
infinitely oscillating between points of minimal contraction and maximal expansion 
has vanishing measure in the set of all possible trajectories it can at the same time 
have non-trivial fractal dimensionality. This phenomenon was first discussed in a 
little bit different terms in the paper by Page 14 . 

III. The cosmological model with the massless scalar field and the "anti- 
slow-roll" perturbation theory 

We begin this paper with the rewriting of the equations (2.5) and (2.6) in terms 
of conformal time 



■n = f dt-^— 
' J a(t) 



because in the conformal time the solutions for the corresponding equations of 
motion in massless limit can be written down explicitly. Thus, we have 

2w'a' 

if" + -2-— + m 2 (pa 2 = (3.1) 

and 

a ' 2 , 47T , 9 47T 9 9 

- T = -l + - T / + - T mV. (3.2) 



3mp 3mp 

Here, "prime" denotes differentiation in respect to rj. At to = Eqs. (3.1) and 
(3.2) turn into 

if H — = (3.3) 

and 

^=- 1+ 3^ (3 ' 4) 

These system of equations (3.3) and (3.4) can be exactly integrated. As in the 

preceding section we shall take as an initial moment of evolution that at which the 
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Universe is in the point of maximal expansion. Then we can get 



a = a y / cos~2r/ 



(3.5) 



— = — tan2?7, 



a 



(3.6) 




3m 2 p 1 
4-7T cos 2?7 ' 



(3.7) 



(here, we have chosen direction "up" for the evolution of scalar field just like in the 
preceding section) and, finally 



These solutions correspond to trajectories which begin in one singularity a = 0, ip = 
±oo go through the point of maximal expansion a = a , <p = <po an d go to other 
singularity a = 0, ip = =Foo (see Ref. 16). In terms of our Fig. 2 one can say that 
all half-plane (a, ip) is filled with the only region 0. 

Now, our purpose is to construct the perturbation theory for the solution of 
Eqs. (3.1) - (3.2) using the solution of massless equations (3.5)-(3.8) as zero-order 
approximation. Apparently, such perturbation theory is predestined for describing 
situations when the kinetic term in the Hamiltonian is much greater then the poten- 
tial term in contrast with the well-known slow-roll approximation (see, for example, 
Ref. 7,19). Thus, can call our scheme "anti-slow-roll approximation". Again we 
choose as a starting moment for our evolution that moment, when the Universe is 
placed in the point of maximal expansion. For calculational simplicity we restrict 
ourselves by a symmetric case <p(0) = 0. It will be more convenient to consider 
instead of the variable a the variable h defined as 




(3.8) 



= 2l 

a ' 



(3.9) 
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When representing the solutions of equations of motion as 



ip = ip (a) + Sip (3.10) 

and 

h = h^+8h, (3.11) 

where ip^ and are given by formulae (3.8) and (3.6) respectively, we can get 
the following equations for Sip and Sh: 

Sip" + 2V/i (0) + 2<5V (0) ' + mV°) V 0) (3.12) 

and 

Air Ait m 2 ,n( a ) 2 n( ) 2 

^.^V^SMt-f-. (3 . 13) 

Substituting into Eqs. (3.12), (3.13) explicit expressions for ip^\ ^°)', and a<°) 
from Eqs. (3.5)-(3.8) we have 



<- // ~ „ <■ ; ~<-, /3top 1 1 9 o 3mi, l + sin2r7 _ , x 

V - 2 tan 277(509' + 2£/iW — — + -m 2 an cos 2m/ — In — - (3.14) 

^ ' Y V 4ir cos2ij 4 ' V 4tt 1 - sin 2t? v 7 



and 



h^Sh = J^rAr^' + ^ in 2 l±^L. (3.15) 

]] 3m 2 p cos2r] * 32 l-sin277 v y 

One can easily see from Eqs. (3.14) and (3.15) that the true small parameter in our 

perturbation theory is m 2 ag. 

Now one can using Eq. (3.15) express Sh through Sip. We have 



y 3mp sin 2?7 32 sin2ry 1 — sin2?7 
Substituting Sh from Eq. (3.16) into Eq. (3.14) we come to the following equation 
for Sip: 

Sip" - 2<V (tan 2rj + ] 

\ cos 2?7 sin 2r] ) 

3m p 99 / cos2t7 , 9l + sin2r7 cos2?7 , l + sin2?7\ 

— J -m 2 ai „ . ' In — — - In — - = 0(3.17) 

4tt V 16 sin 2ti 1 - sin 2r] 4 1 - sin 277/ v ; 
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This equation can be integrated and its solution looks as follows: 



, 3m 2 P m 2 al sin2?7 
1 4tt 16 cos^ 2r] 



l + sin 2 2r7 2 l + sin2?7 , l + sin2?7 2 2 . 

■ In : cos 2rj m : h In (1 — sm 2rj) — In (1 + sin 2?7) 



„, l + sin2f7 /l— sin2?7\ /l + sin2r?\ . . 
+21n21n- —i + 2Li 2 ( '-) - 2Li 2 ( +C\. (3.18) 



2 sin 2rj 1 — sin 2r/ 1 — sin 2i] 

1 + sin 2r] 
1 — sin 2r) 

Here, in Eq. (3.18) Lii(x) is dilogarithm function, which can be defined as a series 

20. 



n— 1 



Li 2 (a;) = 51 

n=l 

The constant of integration C in Eq. (3.18) should be chosen equal to zero to 
provide the satisfaction of the condition 

6h(0) = 0, 

which means that the first-order correction to that can be obtained by substi- 
tution of Eq. (3.18) into Eq. (3.16) should not destroy the initial conditions for our 
evolution, i.e. that the beginning of our consideration coincides with the point of 
maximal expansion. 

The behavior of 5h and 6<p' are presented in Fig. 6a and 6b respectively. We 
see that 6<p' begins since zero value monotonically decreasing and tending to — oo 
at f] — > 7r/4, while Sh beginning since zero value monotonically increases and tends 
to +oo at 7] — ► 7r/4. Both first corrections have opposite signs in respect with the 
corresponding functions in zero approximation. So one can say that inclusion of the 
first corrections Sip' and Sh can describe the transition from the set of trajectories 
without bounces and </?-turns to those, having more complicated structure. It is 
interesting also to depict the plots for the relative values and which are 

Sh 

hi") 



presented in Fig. 6c. The absolute value of is everywhere more than that of 
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It means that the we have bounce earlier than t^-turn and in such a way 
our perturbative scheme describes the transition from the region to the region 
1 described in the preceding section. (It looks quite natural that due to the first 
approximation of our perturbative scheme we can describe only the small part of 
the general set of trajectories, besides the regions disposed to the right from the 
region 1 correspond to the comparatively large values of our parameter m 2 a^ and 
hardly can be described by perturbative methods). 

However, here we should recognize that our perturbative scheme fails in the 
vicinity of singularity. Indeed, it is enough to look at the Fig. 6c to understand 
that -^j reaches the value — 1 at any value of the parameter ma that means that 
we shall have bounce on every trajectory independently of the initial conditions. 
However, we know that at small values of the radius of maximal expansion ao we 
have the trajectories lying in the region 0, i.e. the trajectories without bounces and 
(/3-turns. Thus, we see that we can not, unfortunately, to estimate the parameters 
characterising the change of the regimes of the behavior of our trajectories using our 
perturbative scheme, because this scheme does not work in the vicinity of singularity. 
However, one can try to extract some useful information combining the perturbative 
results, numerical methods and the equations describing the location of possible 
points of maximal expansion and minimal contraction presented in the preceding 
section (see also Ref. 9). 

First of all let us consider the symmetric trajectories, i.e trajectories whose 
points of maximal expansion lies at if = 0. Numerical calculations show that the 
boundary between the regions and 1 lies at 



fnO,Q numerical 2.3. 
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(3.19) 



Now one can try to find this value using zero approximation of our perturbative 
scheme or in other words, using the exact solutions of the equations of massless 
model given by formulas (3.5)-(3.8). Naturally, these equations describe the tra- 
jectories without bounces, however, we can put forward the suggestion that zero- 
approximation becomes invalid and the trajectories can have bounce when these 
trajectories come into the region where points of minimal contraction are localized, 
i.e. where according to Eq. (2.10) 

Substituting into Eq. (3.20) expressions for (p(rj) and a(rj) from Eqs. (3.8) and (3.5) 
we have 



4 V a 

man > , — i , ■ o ■ (3-21) 

V ' 1— sin 2rj 

The value of man separating the regions and 1 can be found from the comparison 
of left-hand-side of non-equality (3.21) with the maximal value of its right-hand-side. 
This value is as follows: 

^0 zero approximation — 1.6. (3.22) 

Now one can try to estimate the same value of the parameter man separating 
the regions and 1 using the first perturbative corrections to a and ip. The logic of 
this estimation will be in a way an opposite to that used in the manipulation with 
zero-order (massless) solutions. Indeed, working in zero-order approximation of the 
perturbation theory we have only the trajectories without bounces. However, we 
know in that in some region of our configuration space we can have bounces in the 
full theory. Thus, fixing the moment when our "zero-order" trajectories come into 
this region we can estimate the value of the parameter man giving the boundary 
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between the regions and 1. In first approximation all the trajectories have a 
bounce at the moment when 5h = —hS°\ but we can again estimate the value (pa at 
the moment of bounce and to check if it satisfies the non-equality (3.20). At some 
small values of mao this non-equality is broken and in such a way we can estimate the 
boundary value of mao- Of course, these calculations are much more complicated 
than those for zero-order approximation and need some help of numerical methods. 
It is interesting that the value of mao parameterizing the boundary between regions 
and 1 is 

maQ first approximation — 

1.8 (3.23) 

and it is closer to the value obtained in the framework of full theory (sec Eq. 
(3.19)) than the value obtained in the framework of zero-order approximation (see 
Eq. (3.22)). Thus one can say that up to some extent our perturbation theory 
complemented by some qualitative considerations describes the transitions between 
different types of trajectories presented in the full theory. We would like to hope 
that the developed here perturbative scheme can be useful in other problems where 
we have opportunity to consider the exactly solvable model which is in some kinship 
with the model under consideration. 
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Captions to Figures 

Fig. 1. On the half-plane a > 0, tp the solid hyperbolic curves are those 
separating Lorentzian region from Euclidean one. Dashed hyperbolic curves sepa- 
rate the possible points of maximal expansion a — 0, ii < (these points can exist 
between two branches of dashed curves) and the possible points of minimal contrac- 
tion a = 0, a > (these points can exist between the dashed and solid hyperbolic 
curves). 

Fig. 2. The structure of the region of the localization of possible points of 
maximal expansion is presented. We consider here the points of maximal expansion 
corresponding to the trajectories which at these points has the direction "up" , i.e. 
tp > 0. 

The region corresponds to the trajectories which after going through the point 
of maximal expansion placed in this region go to the singularity a = 0, tp = +oo. 

The region 1 corresponds to the trajectories which after the going through the 
points of maximal expansion have a point of minimal contraction or "bounce" . 

The region 1' corresponds to the trajectories which after the going through the 
point of maximal expansion have a <^-turn, i.e. the point where tp = after that 
they fall to singularity a = 0, tp = — oo. 

The region 2 corresponds to the trajectories which after the going through the 
point of maximal expansion and subsequent <p - turn have bounce. 

The region 2' corresponds to the trajectories which after going through the 
point of maximal expansion and two subsequent tp - turns fall to singularity and so 
on. 

We did not depicted the regions 1,2,... described in the main text. These 
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regions are included into regions l',2', ... correspondingly. 

Fig. 3. The parts of the trajectories corresponding to the regions presented 
on the preceding Figure are depicted. We trace out for the sake of simplicity of 
the classification only the pieces of these trajectories beginning since the point of 
maximal expansion ignoring their "prehistories" . On the Fig. 3a — 3h we consider 
the trajectories whose points of maximal expansion are placed at the line ip = 0, 
though it is not essential from the general point of view. 

In Fig. 3a is presented the trajectory corresponding to the region falling 
after the point of maximal expansion to the singularity. 

In Fig. 3b the trajectory corresponding to the region 1 is presented. This 
trajectory has a bounce after the going through the point of maximal expansion. 

In Fig. 3c is presented the periodic trajectory having two symmetric points 
of "full stops" a = 0, <p = at the boundary between Lorentzian and Euclidean 
regions and the point of maximal expansion placed just at the boundary between 
the regions 1 and 1' (or, to be precise on the boundary between the regions 1 and 
1. 

In Fig. 3d the trajectory corresponding to the region 1 is presented. This 
trajectory after the point of maximal expansion has ip - turn, almost immediately 
after it - bounce, then the second point of maximal expansion and then it goes to 
singularity. 

In Fig. 3e we have the trajectory corresponding to the region 1' which after 
p - turn falls to infinity. 

In Fig. 3f we see the trajectory corresponding to the region 2 which after the 
p - turn has a bounce already in the lower half-plane. 

In Fig. 3g we see the periodic trajectory whose point of maximal expansion 
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lies on the axes ip = just at the boundary between the regions 2 and 2. 

In Fig. 3h is shown the trajectory corresponding to the region 2' which after 
two tp - turns fall to upper singularity. 

In Fig. 3i is shown the trajectory whose point of maximal expansion is placed 
in the upper branch of the region 1. This trajectory immediately after the bounce 
has (p - turn, then goes through the second point of maximal expansion falls into 
singularity in the lower half-plane. 

Fig. 4. The structure of the subregions of the region 1 is presented. Thin 
regions denoted by letters "b,c,d,e" correspond to the trajectories which have at 
least two bounces. The subregion "a" and other subregions placed between thin 
regions correspond to the trajectories falling to singularity after one bounce. 

Fig. 5. The trajectories corresponding to different subregions of the region 1 
are presented. The trajectories depicted on the Fig. 5a — 5e correspond respec- 
tively to the subregions "a-e" presented in Fig. 4. The trajectory depicted on 
Fig. 5a has only one bounce while the trajectories have at least two bounces. The 
trajectories shown in Fig. 5c and Fig. 5e have a point of full stop (or at least are 
very close to it). 

Fig. 6. In Fig. 6a is presented the dependence of the Sh on the conformal 
time T] while in Fig. 6b the dependence of 6<p' on r\ is presented. In Fig. 6c solid 
line shows the behavior of Sh/h and dashed line shows the behavior of 5<p'/ip'. 
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